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We study vortex lattice structures of a trapped Bose-Einstein condensate in 
a rotating lattice potential by numerically solving the time- dependent Gross- 
Pitaevskii equation. By rotating the lattice potential, we observe the transi- 
tion from the Abrikosov vortex lattice to the pinned lattice. We investigate 
the transition of the vortex lattice structure by changing conditions such as 
angular velocity, intensity, and lattice constant of the rotating lattice poten- 
tial. 

PACS numbers: 03.75.Lm, 03.75.Kk 



Quantized vortices are one of the most characteristic manifestations of 
superfluidity associated with a Bose-Einstein condensate (BEC) in atomic 
gases. By rotating anisotropic trap potentials, several experimental groups 
observ ed fo rmation of triangular Abrikosov lattices of vortices in rotating 
BECs CEEl Microscopic mechanism of the vortex lattice formation has 
been extensively studied both analytically and numerical ly using the G ross- 
Pitaevskii (GP) equation for the condensate wavefunction^'^'^'^'^'^^^^. More 
recently the vortex phase diagrams of a BEC in rotating lattice potentials 
have attracted theoretical attention, since one expects vortex pinning and 
structural phase transition of vortex lattice structures Recently, a 

rotating lattice has been experimentally realized at JILA, makinig use of 
a laser beam passing through a rotating mask Stimulated by the re- 
cent JILA experiment, in this paper, we study vortex lattice structures of a 
BEC in a rotating triangular lattice potential created by blue-detuned laser 
beams. 

We numerically solve the two-dimensional Gross-Pitaevskii (GP) equa- 
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Fig. 1. Density profile (a) and Structure factor profile (b) of BEC in 
the lattice potential. The lattice potential geometry is triangular lattice 
at a/aho = 2.2 and a/aho = 0.65. 



tion in a frame rotationg with anguler velocity Q: 



dt 
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Here the total external potential is given by T4xt(r) 

Mattice(r)5 Lz = ~^^(^^ ~ ^ Component of the anguler mo- 

menum operator, g = A.Ti?i^as/m is the strength of interaction with being 
the g-wa ve scatterin g length, and 7 is the phenomenolocical dissipative pa- 
rameter EESEninSl, The lattice potential created by blue-detuned laser 
beams arranged in the lattice geometry is expressed as 



lattice (r) = ^ Vbexp 
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where r^i,ri2 — '^^i^i + n2a2. We consider the triangular lattice geometry 
with the lattice constant a: ai = a(l,0), a2 = a(— 1/2, a/3/2). Throughout 
this paper, we scale the length and energy by a^o — \/fiJrriuJ^ and Tiuo^). We 
set the dimentionless interaction strength as C = AnagN/hz — 1000, where 
N is the total particle number and hz is a height of the cylinder, and the 
width of the laser beam as cr/aho — 0.65. In our parameter set, the healing 
length is i/a^o — 0.12. 

In our simulations, we first determine the ground-state condensate wave- 
function without rotation by setting = in the GP equation (^J. Fig. 1(a) 
shows the condensate density profile in the lattice potential without rotation. 
We then dynamically evolve Eq. (^J with a fixed angular velocigy ^ and look 
for equilibrium states. In addition to the condensate density profile, we also 
look at the density structure factor defined by 



5(k) = / dm(r)e- 



(3) 
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(a) Fo=0.05 (b) Fb=0.1 (c) ^0=0.3 (d) Vq=0.5 (e) Fo=9.0 



Abrikosov lattice < pinned lattice 

Fig. 2. Density profiles (upper) and Structure factor profiles (lower) of ro- 
tating condensates at a/aho = 2.2, a/aho = 0.65 and ft/uo = 0.70. The 
figure show the vortex lattice structures corresponding to the intensity of 
the lattice potential : (a) Vq = 0.05, (6) Vq = 0.1, (c) Vq = 0.3, (d) Vq = 0.5, 
(e) Vo = 9.0. 



where n(r) = |'0(r)P is condensate density. The structure factor provides us 
information about the periodicity of the condensate density. For triangular 
lattice geometry, there are periodic peaks of regular hexagonal geometry. By 
looking at the position of peaks of the structure factor, we can distinguish 
between the Abrikosov lattice and the pinned vortex lattice. 

Fixing the lattice constant with a/aho — 2.2, we investigated the tran- 
sition of the vortex lattice structure by changing the intensity Vq for various 
angular velocities Q/uq. Fig. 2 shows the density profiles and the structure 
factors for Q/cjq = 0.70. One can see that for weak lattice potentials ((a) 
Vo — 0.05, (6) Vq = 0.1), vortices form the Abrikosov lattice, which is in- 
commensurate with the lattice potential. By slightly increasing the intensity 
((c) Vq = 0.3), vortices start to being partially pinned by the lattice poten- 
tail. For strong lattice potentials ((d) Vq = 0.5, (e) Vq = 9.0), all vortices 
are pinned by the lattice potential. We thus observed the transition of the 
vortex lattice structure from the Abrikosov lattice to the pinned lattice. 

In order to quantify the structural phase transition of the vortex lattice 
structure, we calculate the peak intensity of the structure factor S(K) at 
the lattice pinning point, where Ki = 27r/a(l, 1/V3), K2 = 27r/a(0, 2/^3). 
In addition, we also calculate the lattice potential energy 

£^iattice = J dr^*(r)Mattice(r)^(r). (4) 

In Fig. 3, we plot S(K) and £^iattice against of the intensity Vq. In Fig. 3(a) 
for ^/cjQ = 0.70, we find that £^iattice decreases gradually with vortices be- 
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Fig. 3. Lattice potential energy (•) and peak intensity of the structure factor 
at the lattice pinning point (o) with a/aho = 2.2, a/aho = 0.65 corresponding 
to the angular velocity: (a) Q/cjq = 0.70, (6) Q/cjq = 0.55. Here ^'latticeO and 
S'o(K) is the lattice energy and the peak intensity of the structure factor at 
the lattice pinning point of ground state for the intensity respectivity. 



ing partially pinned, and reaches constant when all vortices are pinned for 
Vo > 0.4. Correspondingly S{K.) increases gradually and finally becomes 
constant when all vortices are pinned. From these results together with 
directly looking at the condensate density profile, we conclude that the 
structural phase transition occurs at Vq = 0.4. Fig. 3(b) shows the re- 
sults for Q/ouq = 0.55. In this case, there is an intermediate domain where 
the Abrikosov lattice and the pinnd lattice coexist. In this domain, the 
vortex lattice structure cannot be categorically determined because of the 
competition between the vortex- vortex interaction and the lattice pinning 
effect. 

In Fig. 4, we map out the phase diagrams of the vortex lattice structures 
against O/cjq and Vq. The lower domain is the Abrikosov lattice domain, 
while the upper domain is the pinned lattice domain. The intermediate 
domain represents coexisting state of the Abrikosov lattice and the pinnd 
lattice. From Fig. 4(a) for a/aho = 2.2, we find that the pinning intensity 
takes the minimum value T^min = 0.4 at Q/ujq = 0.70. Similarly, we find 
^min = 0.7 at ft/uo = 0.83 for a/aho = 2.0 (as shown in Fig. 4(b)), and 
Vmin = 1.0 at f^/cjQ = 0.65 for a/aho = 2.4 (as shown in Fig. 4(c)). 

The minimum values of the pinning lattice intensity Vq for each lattice 
constant can be understood as follows. When vortices form the trianglu- 
lar lattice and undergo rigid rotation with angular velocity the lattice 
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Fig. 4. (1) Phase diagrams for the vortex lattice structures: (a) a/aho = 2.2, 
(6) a/aho = 2.4, (c) a/a^o — 2.0. (2) The angular verocity corresponding the 
minimum pinning intensity: (d). The solid line represents Eq. (5). 



constant is expressed as a function of angular velocity by 

where — h/m is the circulation of a singly-quantized vortex. In Fig. 4(d), 
we plot the lattice constant a/a^o against fi/cjQ giving the minimum pinning 
intensity. We find that it is well fitted by the funcion in Eq. Q . This means 
that the pinning intensity Vq takes minimum value when a^iO) matches 
to the lattice constant a. When this "matching relation" is satisfied, weak 
lattice potential has a stronger pinning effect than vortex-vortex interaction, 
which leads to a sharp transition of the vortex lattice structure. 

In summary, we have studied the transition of the vortex lattice struc- 
tures of Bose-Einstein condensates in a rotating triangular lattice potential. 
We showned that the transition is determined depending on parameters of 
the lattice potential. We also found that the lattice potential has a strong 
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pinning effect when the matching relation is satisfied. For a future study, 
we will discuss the transition of vortex lattice structures for various lattice 
potential geometries. 

We thank S. Konabe and S. Watabe for useful discussions and com- 
ments. We also thank N. Sasa for helpful comments on the numerical simu- 
lations. 
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